THE INVERSE ELECTROMAGNETIC SCATTERING PROBLEM IN 
A PIECEWISE HOMOGENEOUS MEDIUM 

XIAODONG LIU*, BO ZHANGt, AND JIAQING YANG* 

Abstract. This paper is concerned with the problem of scattering of time- harmonic electromag- 
netic waves from an impenetrable obstacle in a piecewise homogeneous medium. The well-posedness 
of the direct problem is established, employing the integral equation method. Inspired by a novel 
idea developed by Hahner [11 j , we prove that the penetrable interface between layers can be uniquely 
determined from a knowledge of the electric far field pattern for incident plane waves. Then, using 
the idea developed by Liu &; Zhang I22| . a new mixed reciprocity relation is obtained and used to 
show that the impenetrable obstacle with its physical property can also be recovered. Note that 
the wave numbers in the corresponding medium may be different and therefore this work can be 
considered as a generalization of the uniqueness result of [20] . 
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1. Introduction. We consider the scattering of time-harmonic electromagnetic 
plane waves with frequency a; > by an impenetrable obstacle which is embedded in 
a piecewise homogeneous medium. For simplicity, and without loss of generality, in 
this paper we restrict ourself to the case where the obstacle is buried in a two-layered 
piecewise homogeneous medium, as shown in Figure 11.11 Note that our method and 
results can be easily extended to the multi-layered case. Precisely, let il2 C M'^ denote 
the impenetrable obstacle which is an open bounded region with a boundary Si 
and let ]R^\ri2 denote the the background medium which is divided by means of a 
closed C^'" (0 < a < 1) surface So into two connected domains ilo and fii. Let V. 
denote the complement of fto, that is, M'^\f7o. We assume that the boundary Si 
of the obstacle 1^2 has a dissection 5*1 = Fi U where Fi and F2 are two disjoint, 
relatively open subsets of Si . 

The electromagnetic properties of the homogeneous medium in VIq are described 
by space independent electric permittivity eg > 0, magnetic permeability /ip > 
and vanishing electric conductivity (Tq = 0. The electromagnetic properties of the 
homogeneous medium in fii are determined by space independent electric permittivity 
ei > 0, magnetic permeability fii > and electric conductivity ci > 0. We define 
the wave number kj in the corresponding medium ilj by fc| = [ej + ioj /uj)^jU!^ with 
Kfcj > 0,3fcj > (j = 0,1). 

For convenience we denote hy T{Sj) and T°'"(S'j) (j = 0,1), < a < 1, the 
spaces of all continuous and uniformly Holder continuous tangential fields equipped 
with the supremum norm and the Holder norm, respectively. Then we introduce the 
normed spaces of tangential fields possessing a surface divergence (see section 6.3 in 
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Fig. 1.1. Scattering in a two-layered background medium 

TdiSj) := {a e T{Sj) : Diva e C{Sj)}, 
r°^"(5,) {a e r"-(5,) : Diva G 

Let T'^{Sj) denote the completion of T{Sj) with respect to the usual L^— norm and 
T^{Sj) denote the tangential fields in T^{Sj) with a surface divergence in L^{Sj) [j = 
0, 1). We denote by v{x) the unit normal vector to a surface at the point x. For a 
closed surface it is directed into the exterior of the surface. For two vectors a, & G C'^ 
we write a ■ b for the scalar product and a x b for the vector product. 

Given four tangential fields Ti,T2 G T°'°'(Sq),Ts, G r°'"(ri) and r4_^r"^"(r2), 
the direct problem consists in finding a solution E,H ^ C^(ilo) D C(rio)j F,G G 
C^(r2i) n C{Th) to the MaxweU equations 

cm\E-ikoH = Q, cm\H + ikoE = Q in f^o, (1-1) 
curlF-iA:iG = 0, curlG + ifcii^ = in l^i, (1.2) 

which satisfies the Silver-Miiller radiation condition 

lim {H xx-rE) = (1.3) 

r— >cxD 

where r = and the limit holds uniformly in all directions x/r, the transmission 
boundary conditions 

V X E - \ev X F = Ti, V X H - Xhv xG = T2 on 5o (1.4) 

with constants A^; and \h given by A^; = eo/(ei + ici/w), \h = \/ /io//ii, and the 
boundary condition 



SS{F) = on Sx 
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with the operator SS depending on the nature of the obstacle $72- Precisely, the 
boundary condition on S'l is understood as: 

vy.F = T-i onTi, (1.5) 

1/ X G- X F) X i/ = T4 onTa. (1.6) 
fci 

with a positive constant A. Note that the case Fa = corresponds to a •perfectly 
conducting obstacle and the case Fi = leads to an impedance boundary condition 
corresponding to an obstacle which is not perfectly conducting but does not allow the 
electromagnetic wave to penetrate deeply into the obstacle. The scattering problems 
with mixed boundary conditions widely occur in practical applications, e.g., in the 
use of electromagnetic waves to detect " hostile" objects where the boundary, or more 
generally a portion of the boundary, is coated with an unknown material in order 
to avoid detection. We refer to [Sj O [6] for the physical relevance and practical 
implication of the electromagnetic scattering by obstacles with a mixed boundary 
condition in a homogeneous medium. 

In the next section, an integral equation method is employed to establish the 
well-posedness of the direct problem. A mixed reciprocity relation will also be proved. 
These results will play an important role in the proof of the uniqueness results in the 
inverse problem. 

The radiation condition (II. 3|) ensures uniqueness of solutions to the exterior 
boundary value problem and leads to an asymptotic behavior of the form 

E{x)^ ——\e^{x) + 0{--)\ , aslxl^oo (1.7) 
fI I \A J 

uniformly in all directions x = a;/|a;|, where the vector field defined on the unit 
sphere S"^ is known as the electric far field pattern. 

We consider the scattering of electromagnetic plane waves 

E\x,d,q) = -^curlcurlge*'=°^-'^ = iUd x q) x de^'^"^''^, 
ko 

H\x,d,q) = cnvlqe''""'-'^ = ikod x ge*''"^-'^ 

where the unit vector d describes the direction of propagation and the constant vector 
q gives the polarization. Then we have Ti = —ly x E^{x, d, q),T2 = —v x H^{x, d, q) on 
5*0, Ta = on Fi and T4 = on F2 in the boundary problem (HHJ-dllll). Throughout 
this paper, we will indicate the dependence of the corresponding scattered field, total 
field and far field pattern on the incident direction d and the polarization q by writing 
E''{-,d,q),H''{-,d,q), E{-,d,q), H{-,d,q) and E°°{-,d,q), H°°{-,d,q), respectively. 

The inverse problem we consider in this paper is, given the wave numbers 
kj [j — 0, 1), the constants and A^ and the electric far field pattern E°°{x,d,q) 
for all observation directions x G S'^, all incident directions d € S'^ and all polarizations 
g e K'^, to determine the interface Sq and the obstacle with its physical property 

As usual in most of the inverse problems, the first question to ask in this context 
is the identifiability, i.e. whether an obstacle can be identified from knowledge of its 
far-field pattern. Mathematically, the identifiability is the uniqueness issue which is of 
theoretical interest and is required in order to proceed to efficient numerical methods 
of solutions. 
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For scattering problems in a homogeneous medium, there has been an extensive 
study in the hterature; see, e.g., uniqueness resuhs for scattering from a perfect 
conductor by Cohon & Kress [5] , for scattering from an impenetrable obstacle with the 
boundary condition (II. 6p by Kress [17] , for scattering from a penetrable obstacle with 
transmission boundary conditions by Hahner [11] and for scattering from a penetrable 
obstacle with conductivity boundary conditions by Hettlich f9^. However, few results 
are available for the case of a piecewise homogeneous background medium. For the 
scalar Helmholtz equation, uniqueness results have been investigated in our recent 
paper [21]. Motivated by Isakov's paper [T^, Kirsch & Kress [T^ gave another proof 
of the unique determination of a penetrable obstacle in a homogeneous medium for 
the scalar transmission problem. Based on their ideas, Hahner proved the unique 
determination of the penetrable obstacle in the electromagnetic scattering using a 
novel method. In this paper, we will use Hahner's method to prove that the interface 
So is uniquely determined from the electric far field patterns in Section |31 Under the 
condition that the wave numbers in the innermost and outermost homogeneous layers 
coincide, that is, ko = ki in the Maxwell equations (|1.1|) - (|1.2|) and 5*0 is known in 
advance, we have proved in [2Q] that the obstacle with its physical property can be 
uniquely determined. A main tool used and established in 1201 is a mixed reciprocity 
relation which is important in the proof of the uniqueness result. In this paper, 
we establish a modified mixed reciprocity relation (cf. Lemma 2.5 in [22 for the 
acoustic scattering). Based on this result, we will prove the unique determination of 
the obstacle and its physical property in Section 21 

2. The direct scattering problem. We first show that the direct scattering 
problem has a unique solution. 

Theorem 2.1. The boundary value problem (I j._?p — ^1.6} admits at most one 
solution. 

Proof. Clearly, it is enough to show that E = H = Q in Hq, F = G = in ili for 
the corresponding homogeneous problem, that is, Ti = T2 = on 5o, Ts = on Fi 
and T4 = on F2. Using Green's vector theorem, we have 



/ u X E ■ Hds — v X E - [{v X H) X v]ds 

J So "'50 

= ^E^H / ly X F ■ [{ly X G) X i/]ds 

Jsn 



XeXh / V X F ■ Gds 

'So 



\e\h i v-X-F ■ Gds + XeXh / (curl F -G-F ■ curl G)dx 

J Si Jq.1 

_Xe^ r x\uxF\^ds + i\E\H I {ki\G\^ -Yi\F\^)dx (2.1) 



where we have used the transmission boundary conditions (jl.4l) in the second equality, 
the Maxwell equations (|l.ip - (|1.2[) and the boundary conditions (|1.5p - (|1.6p in the fifth 
equality. Taking the real part of (j2.ip . we have, on noting that Xe^h ~ fco/fcij 
SRfci > 0, Sfci > and A > that 

fco /■ M ^,2 > 2fcoSftfci3fci f 2 



Therefore, by Rellich's lemma [8], it follows that E = H = in Hq. The transmission 
boundary conditions (jl.4l) and Holmgren's uniqueness theorem jl5| imply that F = 
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G = in Oi, which completes the proof. □ 

Denote by <f>j the fundamental solution of the Helmholtz equation with wave 
number kj {j = 0, 1), which is given by 

<^^{x,y) = ^^^^ x,y&R\x^y. (2.2) 

4:n\x — y\ 

For convenience, let Fq denote the interface Sq. Given two integrable vector fields 
a on So and b on Tj and an integral function tp on Tj, {j — 1,2), we introduce the 
integral operators Mij and Nij , respectively, by 

{A'Uja){x)^2 i^{x) X cm\x{a{y)<i>j{x,y)}ds{y), x G Fj 



{Nija){x) — 2h'{x) X curl curl / zy(y) x a{y)<^j{x,y)ds(y), x £ Ti 

Jso 

for i = 0,1, 2;j = 0, 1, the operators Mij and Nij, respectively, by 

{M,jb){x)^2 i^{x) X cm\^{b{y)<l>i{x,y)}ds{y), x € 
Jr, 

{Nijb){x) — 2v{x) X curl curl / v{y) x b{y)^i{x,y)ds{y), x G F^ 



for i = 0,1, 2,j = 1,2 and the operators Si,j and Ki,j, respectively, by 
{S,,,^)ix) = 2 [ ^i{x,y)ij{y)ds{y), x G F^ 
{K,,,^){x) ^ 2 f ^^l^^[y)ds{y), X G F, 

for i = 0,1,2, i = 1,2. 

Theorem 2.2. The boundary value problem jp — (| J.6'P has a unique solution. 
The solution depends continuously on the boundary data in the sense that the operator 
mapping the given boundary data onto the solution is continuous from T°'"(S'o) x 
T°'"{So) X r°'"(Fi) x T0'«(F2) into CO'"(I^) x CO'"(n^) x C°^''{Th) x C0'"(n7). 

Proof. The uniqueness of solutions follows from Theorem 12.11 We now prove the 
existence of solutions using the integral equation method. We seek a solution in the 
form 



E{x) = ^^cm\ a{y)^o{x,y)ds{y) + XecutIcuiI b{y)<^>o{x,y)ds{y), (2.3) 

<^ So So 



ki 



H{x) = -^curl£;(a;) 



= -^curlcurl / a{y)^o{x, y)ds{y) + -^cm\ / b{y)^o{x,y)ds{y) (2.4) 



for X G f^o and 



F(x) = curl / a(y)"I>i(x, y)(is(?/) + curlcurl / b{y)^i{x,y)ds{y) 

•/So ^ So 
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+curl / c(y)$i(a;,?/)ds(y) + -^curlcurl / i^{y) x {S'^c){y)<Pi{x,y)ds{y) 





= — — curlcurl / a{y)^i{x,y)ds{y) — ikicml / b{y)^i{x,y)ds{y) 



+— curl / d{y)'Pi{x,y)ds{y) 




ki 




ki 



iov x £n\ Si, where a, 6 e T°'"(S'o),c € T°'"(ri),d £ T0'"(r2) and V e C°^°'iT2) 



are five densities to be determined and Sj is the single-layer operator given by 



If the densities a,b e Td{So),c € T°'°'{Ti),d G T°^"{T2) and ^ £ C°'"(r2) satisfy 
the following system of integral equations (j2.7p - (|2.1ip . then we can conclude from 
the integral equations p.7p and ()2.8p that a, 6 € T^'"(S'o). Therefore we choose the 
solution space X := Td(S'o) x Td(S'o) x r°^"(ri) x r0'"(r2) x C°'"(r2) for the system of 
integral equations (|2.7p - (l2.1ip below. The vector field F given by (|2.5p clearly satisfies 
the vector Helmholtz equation and its cartesian components satisfy the Sommerfeld 
radiation condition. Hence, if we insist that divF = in K.-^ \ r2 , then by Theorems 
6.4 and 6.7 in [8] we have that F,G satisfy the Maxwell equations (|1.2p . Since divi^ 
satisfies the scalar Helmholtz equation and the Sommerfeld radiation condition, by 
the uniqueness for the exterior Dirichlet problem (see Theorem 3.7 in [8]) it suffices 
to impose div-F = only on the boundary r2. Then the jump relations together with 
the regularity results of surface potentials imply that E, H, F, G defined in (|2.3p - p.6p 
solve the direct problem provided the densities a, 6, c, d, ip solve 






a + Lia + Mib + Nic + Pid + Qiij} 



2Ti 



on So, 



(2.7) 



Afi/co + A/f fci 



b + L2a + Ahb + N2C + P2d + Q2V' 



2T2 



on Si 



(2.8) 

(2.9) 
(2.10) 
(2.11) 



c + Lga + Msb + N3C + Psd + Qsip 
d + L^a + M^b + N^c + P^d + Q^ip 
iXip + P^d + QsV' 



2T3 




on T2 



on r 



on r. 
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where 
Ni 

L2 
N2 

Q2^ 

L3 



fcl 



Mo,o - AeMo,i 



Ml := Xe{No.o ^ No,i)R, 
Pi XeRSo,2 + iXXEMo,2RSi, 



-2Xev{x) X grad / <^i{x,y)ip{y)ds{y) + iXXERSo,2{i'Tp), 

-{No,o - Nq,i)R, M2 := — —A/0,0 —Mo,i, 



iki 
Xh 



—No.iR+^Mo,iRSf, 
iki ki 



P7 



H 



iX 
"fcl 



Mo: 



-A/0,0 - 
fcl 



—2—^v(x) X curl 
fcl 



v{y)'4){y)^i{x,y)ds{y), 



A/1,1, A/3:=iVi4/?, 
-RSi,2 - iXMi,2RSl, 



N3 ■.^Mi^i + -jNi^iPSl 



P3 
Q31P ■■ 
Li 

Ni 

Pi 

Qiip := 2iXi^{x) X / grad^$i(x, y) x {iy{y) - iy{x)}-)p{y)ds{y) + X^ P 82,2 (i^ip), 



2v{x) X grad / <^i{x,y)-ip{y)ds{y) - iXRSi^2{vip), 

N2.1 - iXRM2,u Mi := fc? A/2,1 - iXRN2,iR, 

N2,iR - iM2.iRSl - iXRM2,i + ^RN2.iPSf, 

fcl 

A/2,2 + iXN2,2PSl ~ iXPS2,2 ~ X^RM2,2RSl + X^PSl, 



Phd 



grad^<I>i(a;,y) • d{y)ds{y), 



Q5^kfS2.2+tXK2,2. 



r2 



Here R, P are defined by Ra a x v and Pa [v x a) x i/, respectively. 
Writing the system of integral equations (|2.7p - (|2.1ip in the matrix form 



fXa 



















(Li 


Ml 


Ni 


Pi 


Qi\ 








(2TA 





Xb 













b 




L2 


A/2 


N2 


P2 


Q2 




b 




2T2 








1 










c 


+ 


L3 


A/3 


N3 


P3 


Q3 




c 




2T3 











1 







d 




Li 


A/4 


Ni 


Pi 


Qi 




d 




2Ti 


Vo 








P5 


ix) 








\0 











Q5J 








\ / 



where A^ = A^; + (A//fco)/fci and A^ = —iXEko ~ iXnki, we can see that the first 
matrix operator has a bounded inverse because of its triangular form and the second 
one is compact by Theorems 6.15 and 6.16 in 8 and the argument in the proof of 
Theorems 6.19 and 9.12 in [8]. Therefore the Riesz-Fredholm theory can be applied 
to prove the existence of a unique solution to the system p.7l) - (|2.11l) in the solution 
space X. To do this, let a,b,c,d and ip he a solution to the homogeneous form of 
([2Jl) - (|2lT|) (i.e., Ti = r2 on 5*0, T3 = on Ti and T4 = on r2). We then 
need to show that a^b = c = d = 'ilj = 0. First, by Theorem 12.11 we conclude that 
E = H = in 0.0 and F = G = in ni. 
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Now, by the jump relations we have 



vxF-^c, -V y. G- = X Sic on Ti, (2.12) 

ki 



X sic 


on 












on 




iki 






^-4, 


on 


T2 



3ld, -z/ X G_ = -^d on T2, (2.13) 



-divF_ = -jy • F_ = -tp on T2 (2.14) 

where G_ denote the limit of F, G on the surface from the interior of Si, 
respectively. Hence, using Green's vector theorem, we derive from (|2.12p - (|2.14p that 

i I \Sic\^ds + i\ I \S2d\^ds + iX ( 
JVi 

= / {v xTZ- curl F + v ■ 'FZ(^iYF}ds 
J Si 

= [ {Icurli^p + IdivFp - [(sRfci)2 - (5fci)2]|Fp - 2mki':ski\F\^}dx. 
Jn2 

Taking the imaginary part of the last equation gives that = on S'lC = on 
Fi and S2d = on This implies that c = on Fi, d = on r2 (see the proof of 
Theorem 3.10 in [8 ). 

On the other hand, define 

E{x) = -mcuA / a{y)^o{x,y)ds{y) r — — !-curlcurl / b{y)^oix,y)ds{y), (2.15) 

J So ^Hko J So 

H{x) = — curli?(x) 
tka 

= — curlcurl / a{y)<$oix,y)ds{y) curl / b{y)<^oix,y)ds{y) (2.16) 

Ko JSo JSo 

for X fl and 

F{x)=cm\ a(y)$o(a;,y)ds(y) + curlcurl / b{y)<i>o{x,y)ds{y), (2.17) 

•^5o So 

G{x) = ^curlF(a;) 

= —curlcurl / a{y)^o{x,y)ds{y) - ikiCvA / h{y)^n{x,y)ds{y) (2.18) 
^'^i JSo Jso 



for a: £ r^Oj where m is a constant given by m := \/\h/\e. Then by the jump 
relations we have 

~ fci 1 ~ 
vxF — jyxF — a, ly x E -\ v x E — a on 6*0, (2.19) 

vxG-vxG = -ikib, -^^v X H + X H = -ikib on Sq. (2.20) 

Thus, E, H, F, G given by (|2.15l) - (|2.18p solves the homogeneous transmission problem 

curlF-iA:iG==0, curlG + ifciF = in f^o, (2.21) 
curl E - ikoH = 0, curl H + ik^E = in f7, (2.22) 

uxF-—vxE^Q, vxG-^^vxH = {) on ^o, (2.23) 
m m\ E 

lim {H XX - \x\E) ^ (2.24) 
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where the hmit holds uniformly in all directions x/lxl. Similar to the argument as 
in the proof of Theorem 12.11 we find on using Green's vector theorem and by the 
definition of m that 



Therefore, by Rellich's lemma (see Theorem 4.17 in fT), it follows that F = G = in 
r^o- The transmission boundary conditions (j2.23p and Holmgren's uniqueness theorem 
[TS] imply that E — H = Q'mVl. Hence, by the relations (|2.19|) and (|2.20p . we conclude 
that a = 6 = on 5o, which completes the proof. □ 

As incident fields besides the electromagnetic plane waves we are also inter- 
ested in the electromagnetic field of an electric dipole with polarization p which is 
given by 

E\x\Zj,p) = — curlcurl(p<I)(a;, Zj)), W{x;zj,p) = cm\{p^j{x, zj)) (2.25) 



located at the point zj e ftj for x ^ Zj [j — 0, 1). For j — 0, 1, denote by E{-; Zj,p) 
and H{-; Zj,p) the corresponding total waves, by £"*(•; Zj,p) and H''{-; Zj,p) the scat- 
tered waves and by E°°{-; Zj,p) and H°°{-; Zj,p) the corresponding far field patterns, 
indicating the dependence on the location Zj and polarization p of the electric dipole. 

Remark 2.3. In the case when the incident field is given by the electric dipole 
E^{x; zo,p), W{x; zo,p) located at zq G ^q, we have E = E^{x; zo,p), Ti = —v x 
E'^{x; zo,p), T2 — ~v X H'^{x; zo,p), — and T4 = in the boundary value 
problem (|1.1|) — (jl.6|) . In the case when the incident field is given by the electric 
dipole E^{x; zi,p), H^{x; zi,p) located at zi G ili, we have E = E{x;zi,p), Ti = 
Xev X E^{x;zi,p), T2 = Xhv x H^{x] zi,p), T3 = —v x E^{x\zi,p) and ^ —v x 
W{x;zi,p) + -^{v X E^{x; zi,p)) x 1/ in the boundary value problem (|l.ip — (11.61) . 
Note that in the last case E°°{-; zi,p) is the electric far field pattern of the field 
E(x; zi,p) = E'^(x; zi,p) for a; G JIq. This is different from our early work [20 , where 
the total field E{x; zi,p) in Hq is decomposed into the scattered field E^{x] zi,p) and 
the incident field E^{x; zi,p), so E°°{-; zi,p) is the electric far field pattern of the 
scattered field E'^{x; zi,p) rather than the total field E{x; zi,p). 

We now establish the following mixed reciprocity relation which is needed for the 
inverse problem. 

Lemma 2.4. (Mixed reciprocity relation) For scattering of plane waves E'^{x, d, q) 
and electric dipole E^{x\z,p) from the obstacle we have 




Taking the real part of the above equation, we have 




So 



inq- E°°{-d;zo,p) = p ■ E" {zo,d,q), zq G fio, 
Airq ■ E°°{-d; zi,p) = Xe^hP ■ F{zi,d, q), zi G ili 

for all incident directions d ^ S'^ and all polarizations p, g G K'^. 

Proof. Arguing similarly as in the proof of Lemma 3.1 in [20j, we have 

Anq- E°"{-d\z,p) 



(2.26) 
(2.27) 
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{v{y) X E'{y- z,p) ■ H{y, d, q) + v{y) x H^y; z,p) • E{y, d, q)}d.s{y) (2.28) 

So 

for all z e Oo U fii and all d e S^, p,q e M^. 

We first consider the case z zq £ ilo- Using a similar argument as in the proof 
of Lemma 3.1 in [20], we have 

p- E''{z,d, q) 

My) X H{y, d, q) ■ E\y- z,p) - v{y) x H\y; z,p) ■ E{y, d, q)}d.s{y) (2.29) 

So 

for all z e fla and all d e S*^, p, g e R'^. Subtract the equation (|2.29p from the 
equation (|2.28p to obtain that 

inq ■ E°^{-d] z,p)-p- E%z, d, q) 

My) X E{y; z,p) ■ H{y, d, q) + v{y) x H{y; z,p) ■ E{y, d, q)}ds{y) 

So 

= AbA// / {i^{y) X F{y; z,p) ■ G{y, d, q) + v{y) x G{y; z,p) ■ F{y, d, q)}ds{y) 
J So 

= ^E>^H / {i^iy) X F{y; z,p) ■ G{y, d, q) + v{y) x G{y; z,p) ■ F{y, d, q)}ds{y) 

JSi 

= 

for all z d 0,0 and all d e S*^, p, g G R'^, where we have used the transmission boundary 
condition (jl.4p in the second equality, Green's vector theorem and Maxwell's equations 
()1.2p in the third equality and the boundary conditions (jl.Sp and (|1.6p in the fourth 
equality. 

We now consider the case z := Zi G fii. Using the transmission boundary condi- 
tion (11.41) and Green's vector theorem, we deduce that 



Anq- E°°{-d;z,p) 

= AijAff / {i^iy) X E'{y; z,p) ■ G{y, d, q) + v{y) x H\y- z,p) ■ F{y, d, q)}ds{y) 
J So 

+XeXh / {i^{y) X F{y;z,p) ■ G{y,d,q) + !^{y) X G{y;z,p) ■ F{y,d,q)}ds{y) 
J So 

= A^Ah / {i^iy) X E\y- z,p) ■ G{y, d, q) + v{y) x H\y- z,p) ■ F{y, d, q)}ds{y) 
J So 

+\e\h I {i^{y) X F{y;z,p) ■ G{y,d,q) +iy{y) X G{y;z,p) ■ F{y,d,q)}ds{y) (2.30) 

JSi 

for all z G ill and all d E , p, q G S."^ . From the Stratton-Chu formula (see Theorem 
6.2 in [8]): 



F{z, d, q) -curl / v{y) x F{y, d, g)$i(z, y) 
J So 

+ curlcurl / vly) x G{y,d,q)^i{z,y)ds(y) 
Jso 

+curl / i^{y) X F{y,d,q)'^>i{z,y) 
J Si 

-— curlcurl / i^{y) x G{y,d,q)<i>i{z,y)ds{y), 
Jsi 
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and with the help of the vector identities 

p • curUcurl^(a(y)$(z,?y)) = a{y) ■ cnT\yCVLT\y{p^{z,y)), (2-31) 
p-curU(a(y)$(z,y)) =a(y)-curlj,(p$(z,?/)), (2.32) 

it foUows that 

Xe>^hP ■ F{z,d, q) 

= -AfiAff / {i^{y) X F(y, d, q) ■ H'{y; z,p) + v{y) x G{y, d, q) ■ E'{y; z,p)]ds{y) 

J So 

+XeXh / {i^iy) X F(y, d, q) ■ H\y- z,p) + v{y) x G{y, d, q) ■ E\y- z,p)}ds{y) 

JSi 

for ah z G rii and all d £ S'^ , p, q £ M.^ . We now subtract the last equation from the 
equation (|2.30|) to obtain 

Anq ■ E°°{-d; zi,p) ~ XeXhP ■ F{zi,d, q) 
= XeXh / My) X [F{y; z,p) + E\y; z,p)] ■ G{y, d, q) 

JSi 

+iy{y) X [Giy; z,p) + H\y- z,p)] ■ F{y, d, q)}ds{y). 
This, together with the boundary conditions 

i^{y) X [F{y; z,p) + E\y; z,p)] i^{y) x F{y, d,q)=0 on Ti 

and 

J^iy) X [G{y;z,p) + W{y;z,p)] - ^[i^iy) x {F{y; z,p) + E'{y; z,p))] x iy{y) 
= t^(y) X G{y,d,q) - -^Hv) x F{y,d,q)] x u{y) =0 on Ts, 

gives the desired result (|2.27p . □ 

3. Unique determination of the interface Sq. In this section we will assume 
that 3fci > 0. Denote by Sq and 5*0 two different interfaces which yield the same elec- 
tric far field patterns for all plane incoming waves, that is, E°°{x, d, q) — E°°{x, d, q) 
for all x, d G 5*^ , q G M^^ . Let and Hq be the bounded and unbounded domains with 

interface 5*0, respectively. We denote by B a large ball containing U SI and by fig 
the connected component of fio H Qq that contains the exterior of i3. If Sq ^ Sq, we 
may assume without loss of generality that there is a point a;* G Sq such that x* G SIq 
and X* G dfl^ (the case with x* G Sq and x* G SIq can be treated similarly). Assume 
that C i? is a ball centered at x* with sufficiently small radius Si > such that 

Bs^ n = 0, as shown in Figure [Ql For fixed 6i we have Bg-^ n SIq + 62i'{x*) C Sle 
and Bsi O Sq — d2i^{x*) C SI for all sufficiently small 82 > 0. 

Lemma 3.1. // the electric far field patterns E'^{x,d,p) for Sq and E°°{x,d,q) 
for Sq coincide for all x,d,q G 5^, then the electric far field patterns also coincide for 
the incoming wave 

(a;) 7^ curlcurl / s{y)^Q{x,y)ds{y), 
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Fig. 3.1. Two different interfaces 

H\x)=cm\ [ s(y)$o(x,?/)ds(y), xeBnile, (3.1) 

■'So 

where s S T^'°'{So) has a compact support in Sq n B^^. 

Proof. By Rellich's Lemma [5], the assumption E°°{x;d,q) — E°°{x;d,q) for 
all x,d e S^, q e implies that E^{z,d,q) = E^{z,d,q) for all z e fie and all 
d ^ S'^, q e R'^. Using the mixed reciprocity relation p.26p . we obtain that the 
electric far field patterns corresponding to both interfaces coincide for incoming waves 
of the form 

1 

E'^(x;z,p) = — curlcurl (p$(a;, z)), 
ko 

W(x;z,p)=cm\{p^Q{x,z)), xeile,x^z. (3.2) 

Furthermore, i?' and i/' defined as in (|3.ip are continuous in B n^e and can be 
uniformly approximated, respectively, by fields of the form 

EsJx) ^ ^cmlcml / s{y)<Po{x + S2iy{x*),y)ds{y), 

Hlix)^cuTl [ s{y)Mx + S2i^ix*),y)dsiy), x eBTwTe (3.3) 

"'So 

by making ^2 sufficiently small. In fact, for a; G (B H ilg) \ Bg^ the differences between 
El^{x) and E'^{x) and between Hg^{x) and H^[x) converge to zero as (52 — >■ since 
Bs^ n 5'o — 52v{x*) C for all sufficiently small 82 > 0, so the kernels are smooth; for 
X G Bg^ the differences converge to zero since E^ and are uniformly continuous in 
i? n fie and Bsi fl fio + S2'^{x*) C fie for all sufficiently small 62 > 0. 

Finally, with the help of a quadrature formula, the integrals in the definition of 
Eg^ (x) and Hg^ (x) can be uniformly approximated on B n fie by sums of waves of the 
form (|3.2p . Note that by the well-posedness of the direct problem we conclude that 
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the electric far field pattern depends continuously on the incoming wave. Based on 
this result we obtain the assertion of the lemma. □ 

We are now ready to prove the main result of this section. 

Lemma 3.2. Let Qfci > 0. // the electric far field patterns E°°{x,d,q) for Sq 
and E^(x,d,q) for Sq coincide for all observation directions x g S^, all incident 
directions d £ and all polarizations g S R'^, then Sq = Sq. 

Proof. We assume that Sq ^ Sq. Then we can choose a point x* as we did at 
the beginning of Section [3l For the same as obtained there, we use a constant 
tangential vector p ^ Q dX x* and a smooth cut-off function p that takes the value one 
close to X* and zero for x with |x — a;*| > ^i/3. Define x„ := x* + di/ni>{x*), n gN. 
Denote by En, Hn, Fn, Gn the solution of the direct transmission problem (|l.ll) - (|1.6p 
for the incoming wave: 



El{x) ^ XEcnilcml / Sn{y)<^oix,y)ds{y), 

JSf, 



where we set 
Then define 

f:{x) 



So 

-ifcoA_ECurl / Sn{y)'^Q{x,y)ds{y), 
J So 



xe Enn,, neN, (3.4) 



s„(x) 



p{x) 



n eN. 



F„ (x) - curlcurl / Sn(y)^i{x,y)ds{y), 
J So 



Gnix) + ikicuA / Sniy)^i{x,y)ds{y), 

'So 



X gQi, n G N 



which satisfies the equations (jl.2p . Since s„ is uniformly bounded in T^(S'o) for all 
n G N, we have that 

T* Jx) -.^i^xEn- Xei^ X f^* = Ab(7Vo,i - NQ^Q)Rsn, (3.5) 
T-l.nix) -.^lyxHn- \ev X = i/coAfiMo^oSn - ikiXnMo^iSn (3.6) 

and DivT* „ are uniformly bounded with respect to the L^(S'o)— norm. Consider the 
integral equations (jOl - dmi) in the space T^{Sq) x T^{Sq) x T°'"(ri) x T0^"(r2) x 
C"'"(r2) with the right-hand sides „(a;) and T2* „(a;) defined as above and 



Tlnix) := ly X F* ^ -Ni^iRsn, xeTi, 

T4*„(x) ■.^iyxG*n- ^{1^ X F„*) XV ^ ikiM2,lSn 

ki 



ki 



RN2ARSn, X£T2 



It can be found that the densities a*,6*,c*,d* and "0* are uniformly bounded with 
respect to the norm of the above space. Using the regularity properties of the integral 
operators in (|2.7p we conclude that the tangential fields a* are even uniformly bounded 
in T^{Sq). Here, the condition S C^'" is important since then we are able to use 
the fact that Mo,o and Mo,i map T^{Sq) continuously into T^{Sq) (see [TBI). Let 
Bsi/2 be the ball centered at x* with radius 6i/2 and define F* :— Sq fl -6^1/2 and 
F 5*0 \ i?5i/2. Then we have that Div Tj* „ is uniformly bounded in L'^(F). Thus, it 
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follows from (|2.8p that Div6* is uniformly bounded in L^{T). These results together 
with the solution representation (|2.5p imply that ly x F* is uniformly bounded with 
respect to the r|(r)— norm. 

By Lemma [Ol and Rellich's lemma £"„ and 77„ coincide with the scattered fields 
from the structure with interface 5*0. Since the incoming fields given by (j3.4p are 
uniformly bounded oii^Oj En and together with their derivatives are uniformly 
bounded in C rio- Hence, by the equation p.Sp and the regularity of „, 

v X F* is also uniformly bounded with respect to the Tj(r*)— norm. 

We now conclude by Lemma 13.31 below that Div {v x G* ) remains uniformly 
bounded in L'^{Sq) for all n e N. This implies that T>\v[Xhv x G„ — x = 
A// Div (i^x G*) + Div [Xhvx (G„ — G* ) — x iJ^] is uniformly bounded in L^(r)— norm 
since s„ is compactly supported in Bg^/^r\So so the kernels in X^vx {Gn — Gn) — v^Hn 
are smooth. Furthermore, The function Div [X^v x Gn — v x H^] = Div (i/ x i?„) is 
also uniformly bounded in L^(r*)— norm since iJ„ and its derivatives are uniformly 
bounded in Bs-^/2- By the vector identity Div {ly x V) = —v ■ curly we conclude that 
the fields 

Div [Xhv X (G„ ~Gl)~vx B\\ 



V)w {-ikxXuv X cvvA \ Sniy)^i{x,y)ds{y) 

J So 

+ifcoA£;i' X curl / Sn{y)^o{x,y)ds{y)] 

J So 

-ifciA//i^ • curlcurl / Sn{y)^i{x,y)ds{y) 
'So 



+ifcoA£;i' • curlcurl / Sn{y)^oix,y)ds{y) 

J So 

remains uniformly bounded in L'^(So). Thus, denoting by $o the fundamental solution 
of the Laplace equation, we have that 



lim • curlcurl / Sn{'y)^oix + ev{x),y)ds{y) 

e->0, £>0 Jg^ 

d f ~ 

= n 1 I / Sr,{v)^n{x + ev{x),y)ds{y) 

is uniformly bounded in L'^{So)- Here we have used the identity curlcurl = —AV + 
VdivV and the fact that fcoA_E — fciA// = fcg(l + fciA|^)/(fciA//) 7^ since /cq, 3fci 
and Xh are positive constants. From this we conclude that 

/„(x) = Divs„(2;) + 2 / Divs„(y) {x,y)ds{y), x e So 

Jso clv{x) 

is uniformly bounded in L'^{Sq). Consequently, Divs„ remains uniformly bounded in 
L^(5o) for all rt e N since the above integral equation is uniquely solvable in L^(S'o) 
and the inverse operator is continuous in L^(S'o). Computing Div s„ and omitting the 
terms that are obviously uniformly bounded in L^(S'o) we have that 

{lyjx) xp) ■ [x-x*) 

\x-x*\^/^ ' ^ ^° 

is uniformly bounded in L'^{Sq). This is a contradiction as can be seen by parameter- 
izing So locally around x*. This ends the proof. □ 
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= 


in VLi 




on 


Sq, 




on 


ri, 




on 


T2 



It remains to prove the following lemma. 

Lemma 3.3. Let S/ci > 0. Assume that Ti n S T^'°'{So) is uniformly bounded 
in rJ(S'o), n is uniformly bounded in T^'°'{Ti) and T4 „ is uniformly hounded in 
T^'"{T2) for alln € N. Then there exists a unique solution F*, G* G C^(r2i) nC(r2i) 
to the Maxwell equations 



with boundary conditions 



V X - ^{v X Fl\ X , 

for each n G N. Furthermore, Div (1/ x G*) is uniformly bounded in L'^{So) for all 
n e N. 

Proof. We first prove the uniqueness result. Let Ti n = on Sq, T^.n = on Fi 
and r4_„ = on Then we just need to prove that F* = G* = in Using 
Green's vector theorem and the above Maxwell equations, we have 

= ki [ {curl Gl-F^ + ihF* ■ 'ffjdx 

= ki[ {Gl ■ cuAF^ + iki\F*\'^}dx ~ [ X\i'xF;fds 

= [ {-^|fclG;p + ^[(3f^A:l)2-(3fcl)2]|F*|2-25^fclSfcl|F*|2}da;- / X\t^ x F*\^ds 
Jill JT2 

Taking the real part of the above equation and noting that Sfci > 0, we get F* = 
G* = in 

To prove the existence, we seek a solution in the form 
F*(x)=curl/ a„(y)$i(x,2/)ds(?/) 

"'So 

+curl / c„(?/)$i(a;, y)ds(y) + -2Curlcurl / v{y) x {Slcn){y)^i{x,y)ds{y) 

+ dn{y)<^i{x,y)ds{y)+iXcm\ v{y) x {Sldn){y)'^i{x,y)ds{y) 

+grad / ^/M{y)^i{x,y)ds{y)+iX v{y)%l;n{y)^i{x,y)ds{y), (3.7) 
Jt2 Jr2 

Gni^) = T^curlF„*(a;) (3.8) 

for X G il\Si. The jump relations together with the regularity of the surface potentials 
imply that F*, G* defined in (|3.7p - p.8p solve the mixed boundary problem in fii 
provided the densities a„, c„, dn, ipn satisfy 

- a„ + Mo,ia„ - -^-iVic„ - -^-Pic?„ - -^-Qii/'n = 2Ti,„ on ^o, (3.9) 

AE AE AE 

Cn + LsOn + NsCn + Pddn + Qsljjn ^ 2T3^n OuFi, (3.10) 

d„ + L4a„ + 7V4C„ + P4d„ + (34?A„ = 2r4^„ on Fa, (3.11) 
iAV^n + Psdn + QsV-n = ouFa. (3.12) 
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This system of integral equations is Fredholm in Y T2(S'o) x T°'"(ri) x T"^"(r2) x 
C"'"(r2). Therefore we just need to prove that the corresponding homogeneous sys- 
tem has a trivial solution. Similar to the argument in the proof of Theorem 12. 2[ we 
can prove that V'n = on c„ = on Fi and dn = on F2. Note that F*, G* are 
also well-defined in ilg a-nd are a radiation solution to the Maxwell equations 

curlF* -i/ciC; =0, curlG,; +ifciF^ = in VLq 

with the homogenous boundary condition 

vxF* = Q on 5*0. 

Thus, by Theorem 6.18 in f8| we obtain that F* = in VIq. By the jump relations, 
we have 

-a„ + A/o.ifln = = a„ + Mo,ian on 5o, 

so a„ = on S'q. 

Since the right-hand side of the system of integral equations (|3.9p - p.l2p is 
bounded in Y , then the solution a, c, d, '0 is also bounded in Y . Thus, using the integral 
representation (I3.7P and by Lemma 1 in [TT], Div {v x G*) = —v ■ curlG* = ikiv ■ F* 
is uniformly bounded in L^(5o) for all n S N. This proves the lemma. □ 

4. Unique determination of the boundary 5*1 and its property. In this 
section we will prove that, given that Sq ~ Sq, the impenetrable obstacle fl2 and its 
physical property ^ can be uniquely determined. Note that in this problem we only 
assume that the wave number fci can be any constant with Q/ci > 0, so the result 
obtained in this section is a generalization of that in [20] . 

Lemma 4.1. For ^2,^12 C fl, let G he the unbounded component o/R'^\(Sl2 U 1^2) 
and let E°°{x\d,q) = E°°{x]d,q) for all x,d e S'^, q ^ with E°°{x;d,q) be- 
ing the electric far field pattern of the scattered field E^{x;d,q) corresponding to 
the obstacle U,2 and the same incident plane wave E''{x]d,q). Let z G Qi H G, 
E^ = E^{x;z,p), = H''{x;z,p) and let F — F{x;z,p) be the unique solution 
of the problem 

cmlE - ikoH = cui\H + ik^E = in ^q, (4.1) 
curli^-ifciG^O curlG + ifcii^ = in f^i, (4.2) 
V X E - \ev X F = Xev X E\ V X H - Xrv x G = Xrv x W on 5o, (4.3) 

V X F = -V X E' on Fi, (4.4) 

V X G - ^{v X F) X V = -V X H' + ^{v x E') x v on F2, (4.5) 
fci ki 

lim {H XX- \x\E) = 0. (4.6) 

\x\^oo 

Let F —F{x; z.p) be the unique solution of the above problem with VL2 replaced by ^2 
and r^i replaced by D,i := fi\(rii). Then 

F{x;z,p) = F{x;z,p), V .t G rJi n G. (4.7) 



Remark 4.2. By Theorem 12.21 the problem (|4.ip - (l4.6p has a unique solution. 
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Proof. By Rellich's Lemma [5], the assumption that E°°{x; d, q) = E°°{x] d, q) for 
ah i, d e S*^, q e imphes that 

V X E{x, d,q) = V X E{x, d, q), v x H{x^ d,q) = v x H{x, d, q) 

for ah a; g 5*0 and ah d e 5^, g G R'^. By Hohrigren's uniqueness theorem (see Lemma 
3.2 in [2i or Theorem 4.1.2.4 in [15 ), it fohows that 

F{z,d,q) F{z,d,q), V d G 5^ q G R^ 

For the electric far field patterns corresponding to the electric dipole we have by the 
mixed reciprocity relation (12.271) that 

E°°{d; z,p) = E°^{d; z,p), WdeS^, qe . 

Rellich's Lemma [5] gives 

v X E{x; z,p) = V X E(x; z,p), v x H{x; z,p) = v x H{x; z,p) 

for all a; G 5o, p G R'^. By Holmgren's uniqueness theorem again it is derived that 

F{x;z,p) = F{x;z,p), V a; G l^i n G, p G R^, 

which is the desired result (|4.7I) . □ 

Making use of Lemma 14.11 and the mixed reciprocity relation (j2.27p we can prove 
the following uniqueness result. 

Theorem 4.3. // there are two obstacles and D,2 which lead to the same 
far field pattern for all observation directions and all incident directions at a fixed 
frequency, i.e., E°°{x,d,q) = E°°{x,d,q) for all x, d G 5^ and all q G M.^ , then 
Si = Si and ^ = SS, that is, the impenetrable obstacle with its physical property are 
uniquely determined. 

Proof. Let G be the unbounded component of \ {0,2^ ^2)- Assume that 
5*1 7^ Si. Then, without loss of generality, we may assume that there exists a point 
zq ^ SiC] ^R'^\r27) • We can choose an /i > such that the sequence 

h , . 

Zj := zq + -v{z^), J = 1,2, ... , 
] 

is contained in G, where 1^(2:0) is the outward normal to Si at zq. Consider the solution 
to the problem (|4.ip - (|4.6p with z being replaced by Zj. By Lemma WA\ it follows that 
E'^{x; Zj,p) — E^{x; Zj,p) for all a; G G and all polarizations p G R^. Since zq has a 

positive distance from fl2, we conclude from the well-posedness of the direct problem 
that there exists a G > such that \^{E''{zo': Zj,p))\ < C uniformly for j > 1 and all 
polarizations p G R'^. On the other hand, by the boundary condition on Si, 

mE'{zo;z,,p))\ = \mE'izQ;Zj,p))\ = \-^{E\zo,z„p))\^^ 

as j ^ 00 for pJ-vlzo). This is a contradiction and therefore Si = Si. 

We now assume that the boundary conditions are different, that is, !^ ^ SS. Since 
the obstacles and the far field patterns are the same, we have F{x,d,q) — F{x,d,q). 
Define 5I2 := ^^2 ~ ^2 and F := F{x,d,q) = F{x,d,q) and consider the case of 



18 



X. LIU, B. ZHANG & J. YANG 



impedance boundary conditions with two different positive constants A 7^ A. Then 
from the boundary conditions 

X G ~ ^{ly X F) X 1^ ^0, vxG-^(vxF)xv^Q ori Si, (4.8) 
ki ki 

we observe that 

{\-\)vxF = Q on Si. 

This, together with the boundary conditions (|4.8[) . implies that vxF^vxG^Q 
on Si. Therefore, by Holmgren's uniqueness theorem [U[TS], F = G = in ili. Using 
Holmgren's uniquness theorem again and with the help of the transmission boundary 
conditions, we conclude that E{x,d,q) = E^{x,d,q) + E''{x,d,q) = in V,q. The 
scattered field E^{x,d,q) tends to zero uniformly at infinity, while the incident plane 
wave E^{x,d,q) has modulus |fco((i x g) x d\ everywhere. Thus the modulus of the 
total field tends to \kQ{d x q) x d\. This leads to a contradiction, giving that A = A. 
The cases with other boundary conditions can be dealt with similarly. □ 

We summarize the main results for the inverse problem in the following theorem. 

Theorem 4.4. Let '<ski > 0. Then the interface Sq and the obstacle ^2 with 
its physical property are uniquely determined by the electric far field patterns 
E°°{x,d,q) for all observation directions x £ 5^, all incident directions d G S^ and 
all polarizations q gM.^ . 

There is a widespread belief that the electric far field pattern for a single incident 
direction d € S^ and a single polarization direction q uniquely determines the general 
obstacle, since the far field data depend on the same number of variables, as does the 
obstacle to be recovered. However, this result remains a challenging open problem [3]. 
Recent progress has been made by Liu, Zhang & Zou [TH] who established uniqueness 
with a single incident wave for a polyhedral obstacle and by Kress [T7] who proved that 
a ball and its boundary condition (for constant impedance A) are uniquely determined 
by the far field pattern for one incident plane wave. In a recent paper llOj . we 
proved uniqueness in determining a small perfectly conducting ball in the inverse 
electromagnetic scattering problem by a finite number of electric far field patterns 
with a single incident direction and polarization. We now extend Kress's result to our 
case that the background is a piecewise homogeneous medium. 

Corollary 4.5. Let Qfci > 0. Assume that the interface Sq and the boundary 
Si are concentric spheres with center at the origin and the impedance X is a posi- 
tive constant. Then they are uniquely determined by the electric far field patterns 
E°°{x,d,q) for all observation directions x (£ S'^ , one fixed incident direction d G S^ 
and all polarizations q (zM.^ . 

Proof. By symmetry, the electric far field pattern for scattering of plane waves 
by the concentric spheres So and 5*1 and the positive constant impedance A satisfies 
E°°{Qx, Qd, Qq) = QE°°{x, d, q) for all x, d € S'^ , q € M.^ and aU rotations Q, that is, 
for all orthogonal transformations with detQ = 1. Hence, knowledge of the electric 
far field pattern for one incident direction implies knowledge of the electric far field 
pattern for all incident directions. The statement now follows from Theorem 14.41 □ 

Karp's theorem in our case as stated in the following corollary is also true. 

Corollary 4.6. Let > 0. Assume that the electric far field patterns 
E°°{x, d, q) satisfies E^{Qx, Qd, Qq) = QE°°{x, d, q) for allx,d€S^,q€ and all 
orthogonal transformations Q with det Q = 1. Then the interface Sq and the boundary 
Si are concentric spheres with center at the origin. 
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Proof. We define 5*0 = Q{So) and 5*1 — Q{Si) for some ortliogonaf transformation 
Q with det Q = 1 . Tlien, by symmetry tlie corresponding electric far field pattern 
E°°{x,d,q) for 5*0, 5*1 satisfies 

E°°{x,d,q) = QE°^{Q-^x,Q-^d,Q-^q) = E^{x,d,q) 

for all X, d G S*^, g G K^. Theorem implies that Sq = Sq and 5*1 = 5*1. This holds 
for all orthogonal transformations Q. Therefore Sq and Si are concentric spheres with 
center at the origin. □ 
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